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I. INTRODUCTION 


Recently there has been growing interest in the issue of parton-hadron duality in QCD 
[l]-[4j]. This is due to both theoretical and experimental progress. On the one hand, we 
witness further improvements in theoretical methods of investigation of low energy properties 
of hadrons, such as the QCD sum rules and Heavy Quark Effective Theory. These methods 
enable one to calculate the hadronic properties directly from QCD, in a model independent 
way. These methods heavily rely on the validity of the parton hadron duality. 

In addition to theoretical advances, one can now compare the predictions of these meth¬ 
ods with the growing amount of the new experimental data. This leads naturally to the 
renewed interest in the fundamentals of these methods, i.e., in the issue of how reliable they 
are. 


One especially important application of the idea of the parton-hadron duality is the QCD 
sum rules 0 (see also Ref. || for review). For the past 20 years QCD sum rules have been 
widely and successfully used to predict masses (the so called two-point QCD sum rules) and 
the coupling constants (three-point QCD sum rules) of different hadrons and their decays. 

The basic procedure in QCD sum rules is the following one: one calculates the physical 
quantity - the polarisation operator of a certain number of currents in two ways. First, we 
calculate the polarisation operator in terms of quarks and gluons, using asymptotic freedom. 
Then, we calculate the same polarisation operator in terms of hadrons, using dispersion 
relations. One then equates the results from the parton model (the so called theoretical part 
of the QCD sum rule) with the sum over the hadron states (the so called phenomenological 
part of the sum rule). Usually, we are interested only in the properties of the lowest lying 
resonance. The contribution of the higher resonances created by the given currents is taken 
into account using the so called continuum model 0. In other words, we approximate 
the hadron spectral density (i.e., the imaginary part of the polarisation operator) by some 
smeared function. The standard approach is to assume that the corresponding smeared 
function is well approximated by the spectral density of the theoretical part of the sum rule. 
The assumed validity of this approach underlies all practical QCD sum rule calculations. It 
is exactly the issue of whether the parton hadron duality holds. 

The issue of the parton hadron duality and the closely related issue of the right model 
for the continuum are not new. The subject was discussed in detail in the early paper by 
Shifman, Vainshtein and Zakharov ||, following the classical discussion by Poggio, Quinn 
and Weinberg || of the parton hadron duality in the case of the e + — e~ annihilation. 

Unfortunately, one cannot rigorously check the parton hadron duality directly in QCD. 
Although many arguments support duality, one cannot tackle the issue without the explicit 
theory of confinement. Consequently, one has to check the hadron parton duality (or, rather, 
its analogues) in simpler models. The existence (nonexistence) of duality in these models is 
a strong argument for (against) the parton hadron duality in real QCD. 

While we are still unable to study the issue of duality directly in QCD, one can learn 
a lot by studying various exactly solvable models, the simplest of which are the quantum 
mechanical potential models. These models have been used to gain insight into the issue of 
the parton hadron duality in the case of two-point QCD sum rules |9].[T5|.[T6| . 

The goal of the present paper is to study the analogue of the parton-hadron duality for 
three-point sum rules in the quantum-mechanical potential models. The issue has become 
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especially relevant recently, dne to the extensive use of three-point QCD sum rules for the 
determination of different parameters of Heavy Quark Physics (see e.g. the review [[H|]). In 
particular, it was found that different QCD sum rules lead to contradictory values of several 
fundamental parameters of Heavy Quark Effective Theory (HQET), like (B\D 2 /(2rriB)\B) 
EH3; these values, in turn, differ from the ones predicted on the basis of the analysis of 
the experimental data |13],[I4|. This requires us to go back and check once again the basic 
assumptions behind the QCD sum rule method. 

Our present study confirms the old results [p|. p~5l , p~6| that state that the duality works 
excellently for the two-point sum rules. However, for the same models where the two-point 
duality works excellently, we shall see that we may encounter serious problems in the study 
of the three-point sum rules. Namely, we shall see that not only local duality, but also 
generalised duality (in the sense defined in Ref. G3> are violated in the situation when the 
two-point sum rules work excellently. 

We shall consider in this note three basic models: 

A) The harmonic oscillator. This is the potential model described by the potential 


V(r) = ur 2 /2. 


( 1 . 1 ) 


B) The linear oscillator. This model is the basis of the realistic potential models; 

V (r) = car. (1.2) 


C) The last model to consider is the linear potential perturbed by the coulombic interaction 
that imitates the effects of the a s corrections in the potential models: 


V (r) = ar — (3/r. 


(1,3) 


The three-point sum rules for the harmonic oscillator were already studied in Refs. 

(The numerical mistake made in Ref. |17| in the three-point sum rule for oscillator (but not, 
of course, for the QCD sum rules for the slope of Isgur-Wise function) was improved in Ref. 
|H|). In Ref. [|T7|, it was shown that the duality may hold for the harmonic oscillator in 
generalised sense only: a one-dimensional integral of the phenomenological spectral density 
is dual to the corresponding integral of the theoretical part of the sum rule. However, later 
investigation |I8| showed that even this duality does not take place, and the three-point 


sum rules for harmonic oscillator do not reproduce the true values. Here, we shall see that 
the harmonic oscillator case is not an exception but the general situation. The duality in 
three-point sum rules breaks down because the coupling signs become alternating. In all 
three models in question, the duality breaks for the sum rules determining ground state 
matrix elements of the following operators: 


0 \ 


a, 


r. 


(1.4) 


(Note the close analogy between the operator ()\ and the operator that determines the 
fundamental parameter of the HQET - the slope of the Isgur-Wise function). 

As a result of the duality breaking, the standard continuum model 0 does not describe 
the true spectral density for the first several resonances. So, the sum rules give answers that 
differ from the right ones (here, in quantum mechanics, we know, of course, the exact values) 
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by 30 — 50%. On the other hand, in the case of the kinetic energy operator 0 2 ~ — d 2 , all 
nondiagonal transitions give positive contributions. The sum rules do work and the duality is 
not broken. Nevertheless, there is a big continuum contribution leading to large uncertainties 
in the predictions. 

This paper is organised as follows. In section 11, we review the notion of duality and 
discuss the duality in the two-point sum rules for the potential models mentioned in the 
Introduction. In section III, we study the duality in the three-point quantum mechanical 
sum rules. The duality fails for the operators O i, O 3 , and we trace the origin of its failure. On 
the other hand, for the operator 0 2 the duality holds, and we investigate the corresponding 
sum rules. In section IV, we give our conclusions and discuss possible implications of our 
results for QCD. The details of the exact solution of the quantum mechanical models at 
hand are given in Appendix A, while the details of the calculations of the theoretical part 
of the sum rules are given in Appendix B. 


II. DUALITY AND TWO-POINT SUM RULES 


Let us recall now in more detail what do we mean by quark hadron duality in the case of 
the QCD sum rules. Consider the function f(q 2 ) at some tensor structure in the polarisation 
operator of two currents in the sum rule. One calculates / by means of the Operator Product 
Expansion in the Euclidean domain of momenta q 2 < 0: 

/ M) ~ a o + a i/q 4 + ■ ■ ■ ■ ( 2 . 1 ) 


Here the coefficient a 0 corresponds to the perturbation theory and the coefficients a^i > 1 
correspond to the matrix elements of the relevant operators over the QCD vacuum. (Possibly, 
the coefficients ao,cq depend logarithmically on q 2 ). The function / is a smooth function of 
its argument q 2 . Its imaginary part, denoted u p (s), is a smooth function as well. Another 
way to calculate / is to express it in terms of hadron properties by means of the dispersion 
relation: 



Im f(s) 

ED 5 )' 


( 2 . 2 ) 


The imaginary part <Jh(s) = Im/(s) involved equals the sum of the delta functions over 
hadronic resonances contributing into the tensor structure in question: 


M s ) = J2M(s~m 2 n ). 

n 


(2.3) 


Now, the parton-hadron duality assumes that, starting at some threshold sq, the integrals 
of the hadronic spectral density ah(s) and the partonic one, a p (s), with exponential weights 
coincide: 

f°<T p (s)e-* /Ml ~ f (2.4) 

J Si Js 1 

The minimal sq satisfying this equation is is called the duality threshold. For QCD sum 
rules to work, si must lie between the masses of the ground state and the first excited state 
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contributing to /. In other words, the hadron spectral density, averaged over some interval 
of s, must be approximately equal to the spectral density calculated in the parton model. 
Of course, this is true for the exact spectral densities. However, our cr p (s) is only a part of 
the exact spectral density, namely the part corresponding to the first several terms in the 
OPE, analytically continued to Minkowsky space. 

Once we know what the duality means in QCD, let us consider its quantum mechanical 
analogue. The quantum mechanical analogue of the polarisation operator of two currents in 
QCD is the time dependent Green function: 

S 0 (0,T) =J2\MV)\ 2 e~ EnT . (2.5) 


Here n runs over the S-states only. This function describes the propagation from the point 
(0,0) to the point (0,T) in the Euclidean time. One can calculate the function Sq in two 
ways. First, we can calculate it as a power series in T for short times. The leading term 
in such a series is an analogue of the parton model contribution in QCD; the higher terms 
imitate the matrix elements of operators in the operator product expansion. Second, one 
can use the explicit formula (|2.5|) and calculate S 0 as a sum over the hadron states. In 


practice, we are interested in the properties of the ground state. So we represent 


*So(0, T) = |'0o(O)pe 


2 o -e 0 t 


+ / dEcr h (E)e 
J E c 


-ET 


( 2 , 6 ) 


where an(E) is the exact spectral density, 

a h (E) = J2^E-E n )\M0)\ 2 - (2.7) 

n=l 

Note that the leading term in the perturbation expansion of the Sq(T) can also be represented 
as the integral of the spectral density cr p (E) times e~ ET . Then the quantum mechanical 
duality means that, after average over some energy interval, 


a h (E) ~ <r p {E), E > E c . 


( 2 . 8 ) 


Equivalently, the integrals of the spectral functions times e~ ET are approximately equal as 
functions of T. The integral of the exact spectral density in the r.li.s. can be also calculated 
as the difference between the exact Green function S 0 (T) and the known exact expression 
for the ground state. Define 


Co(T) 


go(0,r)-|V>o(0)| 2 e- EoT 

|'0o(O)| 2 e-- B ° T 


(2.9) 


and 


C P (T) 


f™dEaJE)e- ET 

i'0o(O)| 2 e' E ° T 


( 2 . 10 ) 


Below the functions C' 0 (T) and C P (T) are called “continuum” functions. If the parton 
hadron duality holds, these two functions must approximate one another for sufficiently 
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small T (corresponding to sufficiently large E). Once duality is established, one can write 
the sum rule to determine the ground state parameters. In order to obtain the sum rule we 
simply rewrite equation (|2.6|) differently: 


roc 

lV’o(0)| 2 e _£ '° T = So(0,T) — / dEa p (E)t 

J E c 


-ET 


( 2 . 11 ) 


Here, for So we use the power expansion, which can be obtained in Perturbation Theory. By 
fitting the r.h.s. of equation (|2.11|) to exponent one is able to determine both ground state 
energy E 0 and its residue 1-00(0)| 2 - The ground state energy can be easily obtained by first 


taking logarithm of (|2.6| ) and then differentiating with respect to time T. Since our method 
is an approximate one, we have to introduce a notion of fiducial region. By fiducial domain 
we mean a window in T where two following conditions hold simultaneously. The first one is 
a control over the power expansion of So- Usually one demands for the last kept term in the 
expansion to saturate less than 30% of the whole expression. This way the upper edge of the 
window is determined. The second condition is a ground state dominance. This requirement 
is needed in order to suppress the relative contribution of the exited states. This condition 
determines a lowest edge of the window. Practically the contribution of all exited states (the 


integral term in ( |2.11|) ) is required to be less than 30%. The fiducial domain corresponds to 
the region where two asymptotics (small T and large T) matches. The important fact is that 
the sum rule (|2.11|) are essentially threshold dependent. Usually we do not know the value 
of the threshold parameter E c except the general point that it should he somewhat below 
the energy of the first exited state (which we do not know too). However, the standard 
philosophy of sum rules is to seek for a region (in E c ) where E c dependence is small. The 
sum rule is then called stable. The variation of the result with E c produces an error, which 
is unavoidable in the sum rule method. Two-point sum rules discussed below happen to be 
very stable with respect to the threshold variation. We depict the sum rules with only one 
optimal value for the threshold parameters, which we define by the best fit to exact known 
results. In all models, the optimal values appear to be very close to the guesses typical in 
practice, e.g., the midpoint between the lowest observed states. 

Let us illustrate the quantum mechanical duality for two-point sum rules for the A, B, C 
models discussed in Introduction [Pj.|Ui|-|TE|| . Although similar discussions are already present 
in the literature, we shall also consider these models for the sake of completeness and as 
simple illustrations of more complicated cases of Section III. 

Consider the harmonic oscillator first. We use the dimensionless units u = 1 and 2 m = 1. 
The left, theoretical part of the sum rule can be represented as a perturbation series: 


^ <0 ’ r > = (Jr) 


3/2 


i _ i T 2 , 2lr 4 - 691 

4 480 120960 


T 6 + • • • 


The corresponding spectral function is 


cr. 


■(E) = 7 -,'TE- 


An 2 


( 2 . 12 ) 


(2.13) 


Consider now the right hand side of the sum rule. The spectral density oy,(U) can 
be represented as a sum of delta functions. For sufficiently high energies the summation 
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can be approximately substituted by the integration [[|.[T7|. Using E n = 3/2 + 2n and 
|^n(0)| 2 = (^) 3/ (2 S )!! , one obtains 


2\ 1 / 2 A 1 / 2 

\ 71 


r / 2 \ U 2 , / 1 \ 3 / 2 

/y B - /2 y 


dE, n 


47T 2 


Ve. 


(2.14) 


The latter is just the theoretical spectral density ( j2.13p . We see that for sufficiently high 
energies the duality holds indeed. 

We can make more detailed estimates, relevant for the sum rules. Define the continuum 
functions Cq(T ) and C£(T) (|2.9|) , ( pTT0| ) . Here we use the exact solutions for the harmonic 

oscillator: U ( l j ar = |; |V , o ar (0)| 2 = ^ . These two functions approximate each other very 

well for T < 1.7, (see Fig. |lj, where the two functions are depicted for the optimum value of 
E c = 2.6). Once duality is established, one can write the sum rule to determine the ground 
state parameters. The fiducial region in this sum rule must, of course, be inside the region 
where the duality holds. Straightforward calculations (see Refs. 0 for details) show that 
this is true indeed, and the resulting values coincide with the exact ones very accurately. 

Consider now the case of the linear oscillator. We use units where a = 1. The l.h.s. of 
the sum rule is given by the asymptotic expansion of the propagator in Perturbation Theory. 
The corresponding Green function was obtained in Ref. [0 (see also Appendix B): 


5 0 l (0,T) = 


(4ttT) 3 / 2 


x _ yy 3 / 2 1 Jl j 13 I ... 
2 12 


(2.15) 


whereas the corresponding parton spectral function is given by eq. ( |2.13| ). Here and below, 
the index L denotes the linear oscillator problem. 

Pass now to the r.h.s. of the sum rule. Proceeding in the same spirit as for the harmonic 
oscillator, the local duality for high energies can be established. We use some exact results 
on the problem, collected in Appendix A. The square of the wave-function |^(0)| 2 equals 
l/(47r). The large n asymptotic behaviour of the energy levels E% = (3/27m) 2 / 3 . The level 
density a%(n) is 


cr 


k 


n) 




n)dn = —. 
47r 


(2.16) 


The energy density is crh(E) = which coincides with the “bare” spectral function 

O- 

The duality for sufficiently high energies established, we turn to the sum rule. The 
appropriate “continuum” functions Cq(T) and C p (T ) of ( fk9l) and ( |2.10| ) are depicted in 
Fig. ^ for the optimal energy threshold E c = 3.4. The fit is perfect and we can plot the sum 
rule for the ground state energy (Fig. |3]). One obtains Eq « (2.35 ± 0.05), while the exact 
value is EQ exact = 2.338. We see the sum rules really do work! 

Our last example of two-point functions, is C model with its “Linear + Coulomb” poten¬ 
tial. The model and its numeric solution are described in Appendix A, while the perturbative 
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expansion of the propagator is presented in Appendix B. Below, numerical solutions of the 
model will be referred to as exact ones. 


^ c(0 ’ T) = WrW l 1 + boV * Tl/2 + hl ^ T 


b n lr 2 


yr 3/2 + ^T 3 + • • •J . (2.17) 


The parameter b 0 = 0.57 is dehned in Appendix A. In QCD, parton spectral density 
acquires corrections due to a s terms in the OPE. In C model, the Coulomb interaction 
imitates the role of these a s terms. The parton spectral function (|2.13|) is modihed: 


a, 


lc 


(E) = — E 1 ' 2 + -^ + —E“ 1//2 . 
47t 2 87t 48 


(2.18) 


Unfortunately, contrary to the previous examples we do not know exact solutions of 
the problem. However, it is natural to believe that the large n asymptotic behaviour of 
the wave-function and the energy spectrum can be obtained in WKB based methods. We 
believe that this way one can confirm that the duality holds in fact. 

Following the procedure described above, the continuum functions Cq c (T) and C’' C (T) 
( |2.9|) and (|2.10| ) were computed and plotted for the optimal threshold E c = 2.9 (Fig. Q). 
Within the window (0.2 < T < 0.9), the duality is valid and we can study the sum rules. Fig. 
D presents the sum rule for the ground state energy. The sum rule result E^ = 1.90 ± 0.05 
matches well the exact one E^ exact = 1.83. 


III. DUALITY IN THREE-POINT SUM RULES 


The goal of the present section is to study the issue of the duality in the three-point sum 
rules. While we have seen in the previous chapter that the duality holds for the two-point 
sum rules, (in quantum mechanics at least), the situation for the three-point sum rules is 
clearly more complicated. Indeed, let us recall the general procedure of the analysis of three- 
point sum rules in QCD ||T§[| . One considers the function f(qf,ql) at the appropriate tensor 
structure of the polarisation operator of three currents. One can calculate this function in 
two ways: using Operator Product Expansion (the theoretical part of the sum rule) and 
saturating by resonances (the phenomenological part of the sum rule). For simplicity we 
shall restrict ourselves here by transitions between the same hadron under the action of 
some current. Then, after the Borel transformation in variables qf, q 2 the phenomenological 
part of the sum rule can be represented as 

f(M 2 , Ml) = ^ 0 2 e- m2(1/M i 2+1/M ^ + J d Sl ds 2 a h (s i, s 2 )er s ^~ s ^ . (3.1) 

Here g is the relevant coupling constant, is the square of the residue of the lowest lying 
resonance created by the current, m is its mass and M 2 , M| are the relevant Borel parame¬ 
ters. Local duality means that the latter integral, taken over some part of the (si, s 2 ) plane 
is well approximated by the corresponding integral of the imaginary part cr p (si,S 2 ) of the 
theoretical part of the sum rule (calculated using Wilson OPE). Even if the local duality 
does not hold, for the three-point sum rules it is possible to have the generalised duality. 
In fact, if there are sign alternating transition, it was argued in Ref. 0 that it may be 


7 









senseless to speak about local duality. The parton model density is likely to be concentrated 
in the narrow area around the diagonal of the (si,s 2 ) plane, while the hadron density is 
spread over the whole plane. It was shown in Ref. |Hj that in this case only “generalised” 
duality makes sense: the partonic spectral density, integrated in the direction orthogonal to 
the diagonal is approximately equal to the hadron spectral density, integrated in the same 
direction: 

/(M 2 ) = gPle~< {1/M2) + r dsa h {s)e ~ s/M2 . (3.2) 

Jso 

Here ah(s) is given by the integral 




(3,3) 


A = (si — s 2 )/2. The parameter so is the continuum threshold (see Fig. |6|). The Borel 
parameters are M 2 = Mf = Mf and we stick to the symmetric point. We can define the 
parton spectral density exactly in the same way: 


/(A/ 2 ) = J 


cr p (s)e S ^ M ds 


(3.4) 


Only for such sum rule there is a hope that the duality (defined in this generalised sense) is 
not be violated. 

Consider now the quantum mechanical analogue of the sum rule (|3.2|) . The analogue of 
the polarisation operator of three currents in quantum mechanics is the function 


<S)(ti,t 2 ) = 


, Ti + t 2 , r, Ti)Oj(r )K(r , n, 0, 0) = 


d 3 rJ2 e ElT2 ^i(®Wi( r )Oi(r)Y^,e £ " Tl ^ n (r)^*(0). 


(3.5) 


Here K(r, Ti,0, 0) is the amplitude of the quark propagation from the point (0,0) to the 
point (r, Ti) in the Euclidean time. At the point (r, T\) the operator () l is inserted. 

The vacuum expectation value (VEV) of the operator Oi is defined as 


(0|6 ( |0) = / d 3 rip'„(r)di(r)ip a {r). 


(3.6) 


The corresponding sum rule can be written as: 

„ roo 

Si(r 1 ,T 2 ) = |^o(0)| 2 (0|0 ? ;|0)e _i?o(ri+T2) + / ds 1 ds 2 cr h (s 1 , s 2 )e~ slT1 ~ S2T2 . (3.7) 

Here, a h (s i, s 2 ) is the exact spectral density. By cr p (si, s 2 ) we denote the theoretical spectral 
function obtained, as for the two-point functions, from the T-expansion of S t . From above 
it is clear that we must study the symmetric sum rules, i.e. T\ = r 2 = T/2 (see Ref. 

for details). We need to compare the cr p (s) integrated with the weight e~ sT , with the 
corresponding integral of (This). Here s — (si + s 2 )/2, and <Jh,p(s ) are the spectral densities 
obtained from the spectral densities cr^ p (si,s 2 ) after the integration over the variable A = 
(si — s 2 )/2 in the same way as in eq. (|3.3|). 

Let us now consider the sum rules and duality for three models considered above and for 
the operators: 



Oi(r) = r 2 /6; 0 2 (r) = — <9 2 ; 0 3 (r) = r. 





A. Harmonic oscillator 


Since the operator O 3 has no analogue in QCD, the sum rules only for two operators 0\ 
and O 2 will be discussed. These sum rules were already investigated in Refs. ||l7l , |l8f . For 
the sake of completeness their analysis is included and extended. 


S\™{T/ 2, T/2) = 
R 2 har (T/2, T/2) = 


32t r 3 / 2 T 1 / 2 
3/2 1 


(4tt) 3 /2 T 5 / 2 
The corresponding parton spectral density is 

1 


1 _ I T 2 ii T 4 _ 692 
3 640 60480 


1 - -T 2 + —T 4 + 
6 


T b + 


og(E) = 
UE) = 


-E~ 1/2 - 
32vr 2 

3/2 / 4 


(4tt ) 3 / 2 ^3^^ 


^ 3/2 _ —E~ 1/2 

6y/7l 


In order to check duality let us, as for the two-point functions, eqs. 
the appropriate continuum functions: 


Sj(T/2,T/2) - |'0o(O)| 2 e-' E ° r (O|Qi|O) 

I V’o (0) 1 2 e~ E ° T (010* 10) 

Ie'c dEcr p (E)e~ ET 


Ci(T) = 

C ■(T) = 

P1V |t/ ; o(0)| 2 e _E ° T (0|Oj|0) 

We shall use here the exact answers for the harmonic oscillator: 


(3.8) 


(3,9) 

|) and ( | 2 . 10 | ), define 


(3.10) 


r„ h "= 3/2; |^S ar (°)| 2 = l/(2ff) 3/2 ; (0|d t |0) k „ = l/3£)„ h "; (0|O 2 |0) h „ = l/2B 0 h “. 

(3.11) 

Let us consider the sum rules for the matrix elements (0|Oi|0)h ar and (0|O 2 |0)har- For 
the operator 0\ both functions C'/ ar and Cp/ r are depicted for the energy threshold E c = 2 
(Fig. 0). An important fact is immediately noticeable. The true continuum is negative and 
cannot be approximated by any positive asymptotics. Hence, no duality persists up to the 
first exited state. Note also that although the continuum contribution to the theoretical 
part of the sum rule is almost negligible (less than 5%), the real contribution of the exited 
states is significant and is about 40% of the ground state. Our choice of E c = 2 was 
motivated by standard guess - it is a midpoint between two observed states E 0 = 3/2 and 
(T 0 + Ei )/2 = 5/2. The displayed picture is not sensitive to the threshold variation and the 
duality is broken for any threshold parameter. 

Let us illustrate how duality breaking becomes fatal for the sum rule. The sum rule is 
obtained by transforming the continuum in equation (|3.7|) to the l.h.s. and then dividing 
the expression by the two-point sum rule (eq. (Q): 
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<0|0i|0> = 


Si(T/2, T/2) — J™ dEa p i(E)e~ ET 


S 0 (T)-f™ 0 dEa p0 (E)e 


-ET 


(3.12) 


Here Si and Sq are obtained by Perturbation Theory. Like in the previous case of two- 
point functions we have to determine a fiducial region. The same conditions of the ground 
state dominance and the control over the power expansions are applied. Of course, for 
the three-point sum rule (|3.12|) we obtain two fiducial domains. One corresponds to the 
numenator (three-point part) and another - to the denominator (two-point part). The final 
window is then obtained in the matching region. In some three-point sum rules discussed 
below continuum contributions are not small. Thus, the standard prescription of the method 
(requirement for the continuum to be less than 30%) may lead to a situation when the window 
almost shrinks to a point. In such cases we increase the bound up to 50%. An important 
notice is that three-point sum rules ( |3.12| ) depend on two (in general independent) threshold 
parameters E® and E c . While E® has to be fitted by the corresponding two-point sum rule, 
E c is a varying parameter of the three-point sum rules. In practice, one usually takes both 
thresholds equal. Below we present some arguments showing that in reality E c is likely to be 
less than its two-point partner E®. Like for the operator 0\ (see above), in all three-point 
functions, which display the duality breaking, we take for the energy threshold E c the value, 
which is somewhat close to midpoint between two lowest observed states. In all these cases, 
the sum rules appear to be almost non-sensitive to the threshold variation and our main 
conclusions on duality violation remain to be valid. 

Fig. |3|(a) shows the sum rule for (0|Qi |0) w with the energy threshold E c = 2.5. Within 
the window (0.5 < T < 1.6) the answer given by the sum rule is about 45% off from the 
exact one ( |3.11| ). Thus, the sum rule leads to completely wrong prediction. In the Ref. 

18| , it was argued that this failure is due to nondiagonal transitions, which are negative 


and numerically large. These transitions are not sufficiently suppressed and they produce a 
strong influence on the sum rule. Sign alternating nature of the exact spectral density <Jh(E) 
leads to the duality breaking at high energies. To illustrate the point, we include in the sum 
rule explicitly a few low lying resonances. Correspondingly, the continuum threshold E c 
rises: 


dEa p i(E)e 


-ET 


N 

E 

k=1 


ale EkT 


+ 


Ie? 


dEa p i(E)e 


-ET 


(3.13) 


Here, k runs over a number N of the first low lying resonances in equation ( |3.5| ). The 
residues are denoted by a\, while the resonances are ordered by their energy levels E k . For 
the case of the harmonic oscillator exact analytic expressions for the energy levels E k and 
the residues of the interest a\ are known |T7|. The energy E k = 3/2 + k. 


CLi. = 


1 4/c + 3 1 (2k + l)\ 

(2tt) 3 /2 6 2^ (fc!) 2 : 


k 


even 


4 = (-!)- 


1 2k + 3 1 (2k + 1)! 


k — odd. 


(2t r) 3 / 2 3 2 2k (k\) 2 

The sign alternating nature of the exact spectral function is clearly observed. 
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We expect the energy threshold E^ to be of order En - the energy of the last explicitly 


taken resonance. Sum rule with continuum of the form (|3.13|) is depicted for N = 3 on Fig. 
H(b) (E^ = 5). The desired plateau is clearly restored. 

The situation with the operator 0 2 strongly differs from the picture described above. 
Despite the fact that the nondiagonal transitions are not vanishing, they are of the same 
sign as the diagonal. Here we use the equation of motion to obtain residues a|: 


= 


Ei. 


( 2 vr ) 3 / 2 


-a 


k ? 


k 


even; a k = 


(-!) 2 a k> k-odd. 


(The factor 3/2 in front of a\ is due to definition of the operator 0\.) 

Thus, the exact spectral function is always positive. This fact will be shown to be crucial 
for the duality to hold. In order to check the duality, the continuum functions C 2 ar and Cpf r , 
eq. (|3.10|) , are depicted (Fig. [|) for the optimal energy threshold E c = 2. Both functions 
match excellently and the duality is established. Consequently, the sum rule for VEV of the 
kinetic energy operator can be investigated (Fig. E). An important remark is in order. In 
the case at hand, continuum dominates in the sum rule (it saturates more than 50%) and 
the window is practically absent. The obtained sum rule displays strong sensitivity to the 
continuum threshold. Thus, such a behaviour of the sum rule is much like the one obtained 
in QCD JT|. Nevertheless, fitting the energy threshold, the exact result ( |3.11| ) can be easily 
reproduced: (0|O 2 |0)har = 0.5 at E c = 2. At this point we disagree with the conclusions 
of Ref. [ 181 on the sum rule failure. In this paper, the continuum threshold was taken the 


same as for the two-point function E c = 2.5. However, in the three-point function at hand, 
E = 2.5 is the energy level of the first nondiagonal transition state. Consequently, for the 
sum rule the value of the energy threshold eventually has to be taken lower. 


B. Linear oscillator 


With the same emphasise on duality, let us investigate the three-point functions for the 
linear potential. The operator O 3 plays now a role of the potential and it is the virial 
partner of the operator 0 2 . Since no exact propagator is known, asymptotic expansions of 
the three-point functions (eq. |3.5| ) are obtained perturbatively. Details of this computations 
are presented in Appendix B: 


5f(T/2,T/2) 

A 2 l (T/2,T/2) 

A 3 l (T/2,T/2) 


1 1 
32vr 3 / 2 T 1 / 2 




/ 247 
\120 



3 1 

167T 3 / 2 T 5 / 2 



fl6y/2- 19 \ 

V 12 ) 



(3.14) 


1 1 
47T 2 T 



y3/2 1 / 35 
8 / V48 



The corresponding parton spectral densities are: 
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^i(E) 

^A E ) 

^s(E) 
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32tt 2 

3 


E ~ l/2 ; 

f . 4 


I67r 3 / 2 \ 3-^/7? 
1 


£ 3/2 


3^, 


4tt 2 


(3.15) 


The general picture with the three-point sum rules for the linear oscillator is very similar 
to the one of the harmonic oscillator. We start from the operator 0 3 and check the duality 
first. Fig. |ll] shows the continuum functions C 3 and C^ 3 (eq. |3.10|) ; E c = 3.3. Again, 
the true continuum is mostly negative and cannot be represented by the asymptotics. A 
best fit would be reached in the “no continuum” approximation. Consider now the sum 
rule for the VEV of the operator. The window for the linear oscillator is moved to the 
left: 0.3 < T < 0.7. Comparing to the exact numerical result (0|O3|0)g Xact = 1.559, one 
(0|O 3 |0) L = 1.1 ±0.1 yield by the sum rule (Fig. [L2](a)) is 35% smaller. We account here for 
the situation, when the positive diagonal transitions almost cancel the negative nondiagonal 
matrix elements. This results in the fact that the large number of resonances must be 
explicitly taken in eq. ( |3.13| ) in order for the sum rule to be saturated. Fig. [T^(b) presents 
the sum rule at A/” = 11 ( The numerical values for the residues a*, are given in the tables [TJ 
and 0. The energy threshold is E c = 7 that lies between third and fourth energy levels. 

Consider now the operator 0\. The sum rule displays the same problem with duality as 
for the harmonic oscillator. No duality persists up to the first exited state. Fig. [L3| shows 
the sum rule for VEV of the operator together with an improved continuum model ( [hi 3| ) . 
The plateau is restored when N = 11 transitions are taken explicitly. 

Let us turn now to the three-point function with 0 2 operator inserted. The corresponding 
hadron spectral function is positive. As it was argued above, positive spectral functions do 
not cause duality breaking. The case at hand confirms this statement. Although the sum rule 
strongly depends on the continuum threshold parameter, it, nevertheless, yields the correct 
value (Fig. [T3|): (0|O 2 |0) L = 0.81 ± 0.01, compared to the exact (0|O 2 |0)^ xact = 0.779. The 
optimal energy threshold parameter is E c = 2.8, which is significantly below the threshold 
parameter corresponding to the two-point sum rule. 

As a common property of the three-point functions with sign changing spectral functions, 
we see that a large amount of resonances must be taken into account explicitly. In other 
words, no duality is valid in the low energy area. This effect is due to the fact , which was 
already mentioned. The window moves to the left, where exited states are not sufficiently 
suppressed. 


C. Linear ± Coulomb model 


Coulomb term added to the potential improves a bit the situation slightly throwing out 
the window to the right. However, the general picture of the duality breaking still persists. 
We now present our results for the three-point sum-rules in “Linear ± Coulomb” potential. 
The three-point functions ( |3.5|) are obtained by perturbation (see Appendix B): 

Si (T/2,T/2) = 327r 3/2 T i/2 x 


12 















( 1+6 “^ rI/2+ ^' 
S^(T/2,T/2) = 


358T - 6 0 1-575T 2 - - 5 $=) T 322 + g -- ^?) % 


(3.16) 


3 1 L + 6 0 4 „ + ^ r 122 + 6 2 1.97T - 6 0 1.38T 2 - ^4_T 3/2 + i 16 ^ 19 ~j r 3^ . 


I67r 3 / 2 T 5 / 2 


30T 


Sj(T/2,T/2) = 

1 1 


4tt 2 T 


^1 + 6 0 1-679T 1/2 + 6gl.477T - 6 0 1.538T 2 - ^ j T 3/2 + 


35 _ 5tt 
48 _ 64 


rj-16 


The corresponding spectral densities are 


<1 (E) 




^s(E) 


1 


32vr 2 

3 


FT 1/2 ; 


I67r 3 / 2 \ 3 ^/¥ 
1 


4 £ 3/2 + + bl2.22E 1/2 - 


3a/? 


3^ 


- 6 n 0.78E- 1/2 


(3-17) 


47r 2 


(l + M.948^- 172 ). 


The behaviour of the sum rule for the kinetic energy operator is similar to ones of the 
harmonic and linear oscillators. Due to the positiveness of the exact spectral function no 
duality breaking is accounted for. Fig. [T5| shows the sum rule at the optimal energy threshold 
E c = 2.4. The exact value (0|O 2 |0)i < l ar f = 0.972 and it matches well the one given by the 
sum rule: (0|O 2 |0) lc = 0.96 ± 0.02. 

Let us now consider the operators O i and 0 3 . Below we plot two graphs representing 


our results on the three-point functions with duality breaking. Fig. [16] shows the usual sum 
rule and the sum rule with N = 5 explicitly taken resonances (E c = 5.5) for VEV of the 
operator 0 3 . Note that this operator is no more a virial partner of the operator 0 2 . Taking 
five transitions explicitly we restore the plateau at the exact level ( 0 |O 3 | 0 )g Xact = 1-401. Fig. 
T7| closes our analysis. It describes the sum rule and the sum rule with N — 3 explicitly 
taken resonances for VEV of the operator 0\, (0|Oi|0)' 


>£act = 0.399. 


IV. CONCLUSION 

Motivated by the recent work on three-point sum rules in QCD, especially connected to 
Heavy Quark Effective Theory, we studied the two- and three-point sum rules in three dif¬ 
ferent nonrelativistic quantum mechanical models with confining potentials. We have seen 
that though in all cases the two-point sum rules work perfectly well, the three-point sum 
rules may fail. Their predictions for the matrix elements of the operators O 13 may differ by 
30 — 50% from the corresponding true values. The reason of the failure is the breakdown 
of duality. The theoretical spectral function is always smooth and positive. On the other 
hand, we have seen by explicit calculation that the “phenomenological” spectral density is 
wildly oscillating and even has a sign changing component due to the sign changing non¬ 
diagonal transitions. We have seen that, though the diagonal transitions between radial 
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excitations have positive transition constants, the nondiagonal transitions always have neg¬ 
ative sign. Consequently, duality does not work for the first several resonances, even if one 
understands the duality in the a “generalised” sense (i.e., even after the integration in the 
direction orthogonal to the diagonal). The averaged “hadron” density strongly differs from 
the “theoretical” one. It seems that the duality starts to work for energies high enough. 
The corresponding threshold lies near the 3d-4th resonance and depends on the model and 
on the operator. However, this is of no practical interest, because one cannot separate the 
leading resonance contribution in order to apply the sum rules. The standard continuum 
model does not work, and the smooth theoretical spectral density strongly differs from real 
one. For the two-point sum rules the whole situation is quite contrary: the duality does 
work in all known examples (section II). 

We have also seen that the absence of the continuum has no relation to the validity of 
the sum rule. In fact, small continuum contribution may arise from the mutual cancellation 
of positive (diagonal) and negative (nondiagonal) transitions. Such a behaviour is displayed 
by the sum rules for the matrix elements of the operators r and r 2 . Moreover, it seems 
that the absence of the continuum contribution is a general feature of the sum rules whose 
right hand side is contaminated by the sign changing transitions. On the other hand, sum 
rules for the matrix element of the operator — d 2 work sufficiently well, despite the fact that 
they are dominated by large continuum. The duality holds thanks to the positiveness of the 
spectral density. This property is very similar to that for the three-point sum rules for the 


transitions under the action of Hamiltonian, where the duality is also not violated [18 


One may conclude that the duality breaking is a general feature of the quantum mechan¬ 
ical analogue of the QCD sum rules in the case of sign changing nondiagonal transitions. 
Furthermore, our study suggests that the duality still holds if all transitions are positive. 
However, the continuum contribution to the sum rule is likely to be large and the result may 
heavily depend on the continuum threshold. 

Our results were obtained for the nonrelativistic quantum mechanical models. It will be 
very interesting to check if our picture of the duality breaking still holds for the relativistic 
analogues of the models A, B, and C. 

Unfortunately, we do not know yet, what may be implications of our results in real QCD. 
Several conclusions, however, can be reached. First, one must be very careful in presence 
of the sign alternating transitions contributing to the polarisation operators. Second, the 
smallness of the continuum contribution is not always of a good omen, and may occur due to 
complicated cancellations in the right hand side of the sum rule. Nevertheless, the situation 
in QCD may be considerably better. Borel transform in QCD may suppress nondiagonal 
transitions stronger than in quantum mechanics, and this may lead to the smallness of the 
contamination. Perhaps, the degree of contamination depends on the matrix element one 
calculates and the form of the sum rule chosen. 

Certainly, further work is needed; in particular, it is desirable to investigate whether the 
complete relativistic calculation improves the situation with duality. Detail examination 
of whether in the three-point QCD and HQET sum rules there are indeed sign changing 
transitions is required. The problem of their relative suppression must be carefully studied. 

Our work was essentially motivated by big discrepancies among values of the matrix 
element (B\D 2 / (2ms)|H) obtained using different sum rules (see Refs. [|TT1 , [T2| ) Consequently, 
we work with simple nonrelativistic analogues for the potential models of B-mesons and study 
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the sum rules for nonrelativistic analogues of the ffQET operator D 2 /(2ms) (and operators 
related to D 2 / (2mg) by the virial theorem). Our results imply that the value for the matrix 
element of the operator D 2 /{2ms) from Ref. [O is underestimated. On the other hand, the 
situation in quantum mechanics seems to be close to that considered in Ref. Jl2|. ffowever, 


the strong dependence on the continuum threshold implies that the accuracy of the relevant 
sum rule may be quite low. Moreover, in the text, we argued that the energy threshold of 
the three-point sum rule should be taken lower (about 20%, according to our experience) 
than that of the corresponding two-point sum rule. The reason is that the contribution 
of the first nondiagonal transition to the polarization operator is suppressed by the factor 
exp[—(£ 1 + Eq)T/2], while the contributions of the first diagonal transition in the three-point 
sum rule as well as of the first excited state in the two-point sum rule are suppressed by 
exp[— E{T\. Consequently, we decided to reexamine the sum rule of Ref. [|T2|. In the latter, 
the energy threshold for the three-point sum rule (co 0 ~ 1 — 1.2 GeV) was taken exactly the 
same as for the two-point sum rule. We have seen above that this treatment of the sum 
rule may lead to a significant overestimation of the matrix element. We investigated the 
leading-order sum rule (eqs. (3.10) and (3.12) of Ref. |l2|) for the kinetic energy operator 
for various three-point thresholds uq ~ 0.8 GeV. The results of our analysis are depicted on 
Fig. [T|. The sum rule yields the following value for the matrix element of the kinetic energy 
operator: (B\D 2 / {2m B )\B) = —0.3 ±0.1 GeV 2 compared to (B\D 2 / (2m B )\B) = —0.6 ±0.1 
GeV 2 for uq ~ 1 GeV. The obtained value is in a good agreement with the ones obtained in 
T3| . [T1|| . (Of course, this result must be considered not as a QCD sum rule prediction, 

[T3| . [T4|1 and QCD 


Ref. 

but rather as an indication that there is no contradiction between Refs, 
sum rule approach). 

After this research was finished we learned that the quantum mechanical duality was 
studied in Ref. [24] for S —> P transitions. However, no violation of duality was found in 


that type of transitions. 
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APPENDIX A: 

I 11 the first part of this Appendix we present some known information about solutions of 
the three-dimensional Schroedinger equation in the case of linear potential. Three operators 
relevant to our study are defined. Their matrix elements are computed numerically and 
represented in tables. Second part of Appendix is devoted to the same analysis for the 
linear ± Coulomb potential. 
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1. Linear potential. 


We look for exact solutions of the three-dimensional Schroedinger equation: 

[- d 2 + ar]V£ = (Al) 

In the dimensionless variables 

E L n =E L n ja 1 E ; x = ra 1 ' 3 . 

equation ( |A 1 |) takes form: 

[-a 2 + = EM. (A 2 ) 

Since the wave-function of the orbitally exited states are vanishing in the origin, we restrict 
our analysis to S-states only. The S-states of the equation (|A2j) are given by the normalized 
Airy functions: 

-0^ = const x Ai(x — E n )/x. (A3) 

Imposing condition of nonsingularity on the wave-function in the origin, we obtain the 
discrete spectrum: 

Ai(— E n ) = 0. 

Table II presents ten first energy levels of the problem. Taking into account the well-known 
fact (Ref. 0) 

Ai(£ -*■ -oo) ~ sin (J ^ 273 + n/4^ , 

we can determine a large n asymptotic behaviour of the energy levels: 

2 


(E n ) ' = m r, n —> oo. 

o 

An important and very special property of the Airy functions is 


l ^( 0)| 2 = ^, 


47T 


(A4) 


(AS) 


and does not depend on n. However, 


ri(0) = (-1)’ 


>/47T 


Below, we present the exact numerical results for the matrix elements of the following 
operators: 

()\ = x 2 /6; 0 2 = -d 2 \ 0 3 = x. 

The matrix elements are defined: 

M l k 3 = (i\O k \j) = J d 3 x^*(x)d k (x)ijjf(x) (A 6 ) 

Note that by the equation of motion 

M? = Ei8 ij - Ml 3 (A7) 

The lowest quarters of the matrices M\ and M 3 age given by tables [TJ and |ITIj respectively. 
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2. Linear + Coulomb potential 


[-a 2 + v>M° = 


(AS) 


The potential 1/(r) is taken of the form: 


V (r) 


4 

= or- 

3 


a s (r). 

r 


ot s {r) 


2 vr 

9 ln(<5 + 7 /r) ’ 


(A9) 


Here, a s (r) is a running coupling constant and it reflects the asymptotic freedom of the strong 
interaction. In order to pass to dimensionless variables the following change of variables is 
performed: 


x = (2 ma') 1 // 3 r; E lc = E^[2m/ a 2 ) 1//3 . 


In the new variables equation (|AS|) reads: 


[-d‘ 2 + x-b(x)/x\^ = E^; 


(A 10 ) 


b{x) = 


87r(4m 2 /a) 1//3 


27\n(8 + 7(2ma) 1 / 3 /^) 


In Appendix B we present a perturbative determination of the propagator for the given 
problem. Unfortunately, this procedure cannot be performed (at least analytically) if the 
running logarithm is present. Since we do not wish to introduce any uncertainty due to the 
logarithm factorisation, we fix the running coupling constant a s at some value. A special 
choice of the parameters is not important for our analysis, but to be concrete we choose 
some quasi-realistic ones |22f: the string tension a = 0.14 GeV 2 ; 8=2; 7 = 1.87 GeV -1 ; 
for the constituent mass we take m = 0.35 GeV. When related problems are investigated 
within QCD the running coupling constant a s is usually of order 0.3. In our study a s is set 
to be 0.28, which corresponds to some fixed point in space, namely xq = 0.088. In the text, 
a parameter b 0 is used and it is defined as b 0 = b(x 0 ) (|A 10 |) . 

I 11 order to solve equation ( |A10| ), we consider the Coulomb potential as a perturba¬ 
tion. Solution of the nonperturbed problem was described above. The Hamilton ( |A10| ) is 
diagonalized in the basis of the wave-functions (PUT): 


n l 




(All) 


k =1 


In order to obtain the low energy spectrum it is sufficient to take into account only a 
few levels (here N], denotes the number of levels). Table presents eigenvalues computed 
by the numerical diagonalization. 

Products of the wave-functions in the origin are 


Cnm = 


1 N l 

lc *(0 ) = t-EM 

k,k' 


(-1) 


k+k 1 


(A12) 
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Coefficients C^ m appear to be slowly changing monotonic functions of the number of 
levels N l taken into account in equations ( All ) and ( A12|) . In order to extract these 
coefficients N L = 60 was taken and then the results were extrapolated. To illustrate the 
point we plot C^q as a function of N Jj (Fig. [19|). This function slowly approaches (as a power 
law) to its limiting value Cqq = 1.83. This procedure may lead to some small numerical 
underestimation, which we estimate not to exceed a five percent level. Table 0 is devoted 
to the coefficients C] c m . 

As for the linear potential, we present tables of numerically calculated exact matrix 
elements (tables [VJ. |V11| . and 


APPENDIX B: 


In this Appendix, we present a perturbative derivation of propagators and matrix ele¬ 
ments of interest. The potential V(r) = ar — (5/r is considered as a perturbation to the free 
particle motion. To study “pure” linear potential j3 must be set to zero. 

The “free” propagator is given by 


G 0 (r,T,r',t') 


Nr^r exp 


(r - r') 2 
4(T - fj 


(Bl) 


Note that in our units 2 m = 1. This factor can be always restored in final expressions. The 
exact propagator G ex can be expressed in a formal perturbation series: 


G ex (r, T, r', t') = G 0 (r, T, r\ t') - J d 3 s £ drG 0 (r, T, s, r)V(s)G 0 (s, r, r', t') + 


d 3 s / dr d 3 s' / dr , G 0 (r, T, s, t)V(s)Gq(s, r, s', t')V(s')Go(s', t', r', t') + 


(B2) 


In the following analysis we retain only terms up to the second order in perturbation. All 
terms of the order a 3 , (3 3 ,a/3 2 ,a 2 /3 and higher will be systematically omitted. We denote: 


G ex = G 0 + aG[ + a 2 G L 2 + /3G% + f3 2 G% + ol(3(G\ c + GfO; 

G\= (Go|s|G 0 ); G L 2 = (Gf|s|G 0 ); Gf = (G 0 |l/s|G 0 ); 

G%= (Gf |l/s|G 0 ); G 2 c = (Gf |l/s|G 0 ); G^ L = (Gf |s|G 0 ). (B3) 


Here, brackets denote the integration over s and t defined in eq. (|B^). Unfortunately, these 
integrations cannot be done in closed form. However, the time integration can be performed. 
This is done with the aid of the following integrals (Ref. |23fl): 



1 

x 2 y 2 

V^x + y . 

{x + y) 2 ' 

[{t - r)r] 3 /2 exp 

( t-r ) r 

,0/0 CXp 

iw /z xy 
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1 

x 2 y 2 

a/7T 

(x + yf 

{t - r) 3 / 2 r 1 / 2 eXP 

(t-r) r 

■^ exp 
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After the time integration is removed, the angular integration can be easily performed 
as well. The remaining integrals over the radial components are of a nonreducible form, 
and we cannot proceed any further. Finally we arrive at the following expressions for the 
propagator components: 

1 r 


G\{r, T, 0, t') = 
G^ir, T, 0, t') = 


(47t) 3 / 2 2(T — t') 1/2 

1 1 / 


fOO x 2 r 2 r 2 

dxe 4 ( T - 4 ') + e 4 ( T - f ') 


dss 4 


32tt 3 / 2 (T - t') 1/2 Jo r J i 

(r+s(l+io))^ (|i—s|+ii>s)^ POO 


dw x 


(r+s(l+io))' 4 (|r — s|+ids)^ POO / ( r-\-s(l+wy ))^ (|i—s|+inj/s)^ 

e W^T~ - e 4 {T-t') _)_ J (Jyy I e _ e mt=F) 


Gf(r,T,0,i') = 

G?(r,T,0,t') = 


r°° dx 

r 2 (l+2x) 2 


1 0 X 

e 4 (T-to 

i 

1 

I 


X 


(47t) 3 / 2 2(T — d) 1 / 2 

r°° ds r°° dx r 1+2x 
I o r Jo x J \i-x\+x 


dw 


e 4(T-t') _ g 4(T-t') 


G^(r,T,0,O = 


r°° dss 2 r 00 


(47t) 3 / 2 4(T — t') 1 / 2 Jo r J i 


dw 


G^(r,T,0,t') = 


o 

— sl+tuxs)^ \ 

(r+s(l+io))^ (|r —s|+iys)^ 

— e 

4 (T-t') 

+ e 4(T-t') _ g 4(T-t') 

00 dss 2 

POO 

/ dw 

/■°° dx 

— x 

r 

Jl J 

0 X 


/ (r+s(l+iu(l+2a;)))^ (|r— s|+it>s(l+2a;))^ 

(1 + 2x) e 4 ^ zf > - e 


(r+s(l+it;(| 1 — :e|+:e)))^ (|r — s|+ius(| 1 — x\-\-x))^ 


+ (|1 —x|+x) e 4 (t3ft 


4 (T-t') 


In order to obtain equations (|2.15|) and ( p . 1 7|) , the limit r —> 0 is taken. At this limit the 
propagator components ( [BBj ) are completely calculated analytically and results of Ref. ||2T[ 
are recovered. 

Having in our disposal the exact propagator ([B3]), we can proceed in computing the 
matrix element of the operators Oj(r). We are interested only in the following amplitudes: 
the free propagation from the point (0,0), in the point (r, T/2), where the operator 0*(r) is 
inserted, and then there is the free motion to the point (0,T). 

(G ex (0, T ,r,T/2)|4(r)|G ex (r,T/2,0,0)) = 

(Go|Oj|Go) + 2a(Gi |0,;|Go) + a 2 (G'f|Oj|Gf) + 2o?^G\ |Oj|Go) + 2(5{G C V |Oj|(?o) + 

/ 3 2 (Gf |Oi|Gf) + 2/3 2 (G^ , |O J :|G 0 ) + 2af3(G°\O i \G L l ) + 2 a(3(G L 2 c + G^ L |d:|G 0 ). (B4) 
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In the equation (|B4|) only three-dimensional r integration is denoted by the brackets. The 
operator 0*(r) is one of the following operators: 

Oiir) = r 2 /6; 0 2 (r) = -<9 2 = ^ ; 0 3 (r) = r. 

For any matrix element in (|B4] ) r enters as a polynomial times Gaussian. Thus, the r 
integration can be easily done. All remaining integrals have a fractional form and can be 
computed analytically or numerically. Final results of these calculations are formulated in 
the text ( 13.14| , |3.16|) . 
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TABLE I. 
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2.223 

4 

-0.004 

0.018 

-0.116 

2.223 

5.610 


TABLE II. 


(i\03\j)L 

0 

1 

2 

3 

4 

0 

1.559 

0.653 

-0.197 

0.101 

-0.064 

1 

0.653 

2.725 

0.974 

-0.275 

0.134 

im 

-0.197 

0.974 

3.680 

1.248 

-0.341 


0.101 

-0.275 

1.248 

4.524 

1.493 

4 

-0.064 

0.134 

-0.341 

1.493 

5.296 


TABLE III. 


El 

E\ c 

El 

El 

E\ c 

El 

El 

E\ c 

El 

El 

1.828 

3.745 

5.245 

6.551 

7.735 

8.833 

9.865 

10.845 

11.783 

12.684 


TABLE IV. 


22 











































































































































S~l\c 

nm 

0 

1 

2 

3 

4 

0 

1.83 

-1.66 

-1.58 

1.53 

1.47 

1 

-1.66 

1.50 

1.44 

-1.39 

-1.35 

IB 

-1.58 

1.44 

1.37 

-1.34 

-1.31 

IB 

1.53 

-1.39 

-1.34 

1.32 

1.28 

4 

1.47 

-1.35 

-1.31 

1.28 

1.25 


TABLE V. 


<*|OiU)ic 

0 

1 

2 

3 

4 

0 

0.399 

0.367 

0.042 

-0.013 

-0.006 

1 

0.367 

1.363 

-0.874 

0.079 

0.021 

2 

0.042 

-0.874 

2.569 

1.471 
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TABLE VI. 
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FIGURES 



FIG. 1. Harmonic oscillator. Duality for the two-point function. Continuum functions Co iar 
and C p are plotted vs. T. The energy threshold E c = 2.6. 
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FIG. 2. Linear oscillator. Duality for the two-point function. Continuum functions Cq and C p 
are plotted vs. T. The energy threshold E c = 3.4. 




FIG. 3. Linear oscillator. Sum rule for the ground state energy. The dashed line corresponds 
to the exact value UQ exact = 2.338. The energy threshold E c = 3.4. 
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FIG. 4. Linear + Coulomb potential. Duality for the two-point function. Continuum functions 
C l 0 c and C l £ are plotted vs. T. The energy threshold E c = 2.9. 



FIG. 5. Linear + Coulomb potential. Sum rule for the ground state energy. The dashed line 
corresponds to the exact value inexact = 1-828. The energy threshold E c = 2.9. 



the generalized duality. 
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FIG. 7. Harmonic oscillator. Duality for the three-point function. Insertion of the operator 
()\. Continuum functions C Y } iar and Cpf are plotted vs. T. The energy threshold E c = 2. 
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FIG. 8. Harmonic oscillator, (a) - three-point sum rules for VEV of the operator O The 
energy threshold E c = 2. (b) - three-point sum rules for VEV of the operator 0\ with N=3 
explicitly taken resonances. The energy threshold E^ = 5. The dashed line corresponds to the 
exact value (0|Oi|0)g a .[ ct = 1/SEq 31 . 
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FIG. 9. Harmonic oscillator. Duality for the three-point function. Insertion of the operator 
O 2 . Continuum functions C and are plotted vs. T. The energy threshold E c = 2. 
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FIG. 10. Harmonic oscillator. Three-point sum rules for VEV of the operator O 2 . The dashed 
line corresponds to the exact value ( 0 |O 2 | 0 )e xact = l/2F'o ar . The energy threshold E c = 2.. 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 

FIG. 11. Linear oscillator. Duality for the three-point function. Insertion of the operator O 3 . 
Continuum functions and Cp 3 are plotted vs. T. The energy threshold E c = 3.3. 
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FIG. 12. Linear oscillator, (a) - three-point sum rules for YEV of the operator 0 3 . The 
energy threshold E c = 3.3. (b) - three-point sum rules for VEV of the operator O 3 with N=ll 
explicitly taken resonances. The energy threshold E^ = 7. The dashed line corresponds to the 
exact value ( 0 |O 3 | 0 )g xact = 1.559. 
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FIG. 13. Linear oscillator. Three-point sum rules for YEV of the operator 0\ and the sum 
rule with N=ll explicitly taken resonances. The dashed line corresponds to the exact value 
(0|Oi|0)g Xact = 0.486. The energy thresholds E c = 3.3, E^ = 7. 
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FIG. 14. Linear oscillator. Three-point sum rules for YEV of the operator O 2 ■ The dashed line 
corresponds to the exact value (0|O2|0)g Xact = 0.779. The energy threshold E c = 2.8. 
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FIG. 15. Linear + Coulomb model. Three-point sum rules for VEV of the operator O 2 • The 
dashed line corresponds to the exact value (0|O2|0)e Xact = 0.972. The energy threshold E c = 2.4. 
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FIG. 16. Linear + Coulomb model. Three-point sum rules for VEV of the operator O 3 and the 
sum rule with N=5 explicitly taken resonances. The dashed line corresponds to the exact value 
( 0 | 03 | 0 }g Xact = 1.401. The energy thresholds E c = 2.9, E^ = 5.5. 
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FIG. 17. Linear + Coulomb model. Three-point sum rules for VEV of the operator 0\ and the 
sum rule with N=5 explicitly taken resonances. The dashed line corresponds to the exact value 
(0|Oi|0)g Xact = 0.399. The energy thresholds E c = 2.9, E£ = 5.5. 
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FIG. 18. The Ball and Braun sum rule to leading-order as a function of the Borel parameter t 
for different values of the continuum thresholds: a) loq = 1 GeV, u\ = 0.7 GeV; b) ujq = 1.2 GeV, 
oj\ = 0.85 GeV ; c) loq = 1 GeV, ui\ = 0.8 GeV. The dashed line indicates the working region. 



FIG. 19. Cqo as a function of N^. 
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